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Complete Feynman diagram automatic computation systems are now coming of age after 
many years of development. They are made available to the high energy physics com- 
munity through user-friendly interfaces. Theorists and experimentalists can benefit from 
these powerful packages for speeding up time consuming calculations and for preparing 
event generators. The general architecture of these packages is presented and the cur- 
rent development of the one-loop diagrams extension is discussed. A rapid description 
of the prominent packages and tools is then proposed. Finally, the necessity for defining 
a standardization scheme is heavily stressed for the benefit of developers and users. 



1. Introduction 

Since its introduction by R.P Feynman0 in 1949, the diagrammatic technique of 
computing matrix element and consequently, all physics quantities in high energy 
physics (HEP) has been extensively used and has proved to be the most simple, 
intuitive and general method to compute even the most complex processes. For its 
conciseness and its pictorial approach, this method has spread over other research 
field like atomic, nuclear and solid state physics. 

Once the Lagrangian of the theory is selected, a process (defined by the initial 
and final state particles and the order of the calculation) is decomposed into a set of 
sub-processes represented by a diagram (fig.|l|). The real power of the Feynman ap- 
proach lies, through the use of definite rules, on a quasi-mechanical transformation 
of each diagram in an algebraic expression representing its quantitative contribution 
to the process. Finally the total cross-section is obtained by the integration over 
the phase space of the square of the sum of the contribution of all individual graphs. 
Although this approach is straight forward in its principles, the actual computation 
of all but the most simple processes is a quite lengthy task, prone to errors and 
mistakes. In order to avoid these pitfalls, calculus were used to be performed inde= 
pendently by several theorists until the results did fully agree. In 1967 M. Veltmaro 
wrote the first computer language to perform algebraic computation of the trace 
expressions resulting from the Feynman diagrams method. SchoonschipS was the 
first step towards an automatization of these quite involved computations. For 
solving this same HEP computation problems, several other developments started 
at this time: i?edwce§, Mocsj/mcll. Nowadays general purpose computer algebra 
languages have been put on the marketplace. They provide valuable tools to per- 
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Fig. 1. Decomposition of a process in Feynman Fig. 2. Main modules of a typical package 
diagrams 

form many physics calculations^: i?educei, Map/ei, MathematicM. However, t 
original Schoonschip and its more modern implementation, the Form languagel 
(J. Vermaseren) are still prefered for these specialized computations. The tech- 
nique developed by R.P. Feynman is well structured and quite suited to a complete 
automatization. However the idea of building a complete package performing the 
computation of any given process only from general principles appeared to many 
as an enormous task and a never ending enterprise. But several groups out-passing 
these arguments decided to launch such projects. This workshop series has played 
a positive role in the motivation and the building of collaboration between these 
groups. Today, several packages can be used by experimentalists to create (yet with 
some limitations) their own event generator in an almost automatic way. Tools and 
libraries have also been created to facilitate complex calculations. 2. Motivations 



These packages which could have remained at the level of pedagogical toys, have 
attracted a lot of interest because modern HEP research requires the calculation of 
more and more complex processes. 

2.1. Larger center of mass energy 

The center of mass energy of the colliders going up (LEP-II, HERA, LHC and NLC), 
new heavy particle thresholds (W, Z , top, Higgs) are crossed increasing therefore 
the number of diagrams and the level of complexity of each diagram computation. 

tp. Nason's paper in these proceedings 
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2.2. Better accuracy 

High precision experiments require high precision theoretical predictions. e + e~ 
colliders experiments belong to this category and channels bearing small contribu- 
tions must be included in the computation to .cope with the high statistic data. 
For example, the t-channel in the e + e~ — > where intringing event accumu- 

lation for large M 77 was observed in L3. —The contribution of the non-resonant 
t-channels in the W production at LEP-I1LL3, often neglected, must be computed 
below and well above the W-pair threshold. Higher precision also means includ- 
ing loops diagrams as for the luminosity measurement at LEP-I which is limited 
by the theoretical uncertainty of the 2-loop calculation of the Bhabha scattering 
(Ath ~ 0.25%, A exp w 0.07—0.16%). The g-2 experiments requires the computa- 
tion of 4-loops contributionsEj to barely match the experimental precision. Loop 
calculations are more involved than tree level ones as the n-dimensional regular- 
ization technique must be applied to treat the singularities. Furthermore, better 
accuracy imposes often to take into account fermion mass effects, radiative correc- 
tions and polarization effects. 

2.3. More Processes 

Looking for small effects revealing possible deviations from the standard model leads 
to the precise computation of many more processes: the "background processes", of 
no great theoretical importance but whose contribution could dilute or completely 
wash out the expected signal. Furthermore in the so-called supersymmetric models, 
the final state particles are often unstable (q — > q + 7, x + — > x° + W + — ► X°99) 
and therefore the actual processes are rather complex (2^> 4,6 or more particles). 
Due to the particle inflation, more propagators are present. For example, e + e~ — > 
e + e~x°X° includes 704 diagrams not taking into account the x° decay. Moreover, 
there are many possible models and each process may need to be computed for each 
model. 

2.4. Perturbative QCD 

Next to leading order QCD computations are highly needed for a precise comparison 
with the experiments, but the complexity of loop calculation in non-abelian gauge 
theories bridles this expectation. New methods beyond the conventional Feynman 
diagrammatic techniques are being developed including spinor hclicity methods, 
supersymmetry Ward identities or string-based techniquescj. The inclusion of these 
new approaches in automatic computation packages is the next logical step. 

In summary, each diagram computation gets more complex (polarization, mass 
effects, loops), each process gets more diagrams (loops, non-standard models, non- 
resonant channels), for each model, more processes (background processes,..) must 
be computed and finally more models have to be investigated. Therefore the need 
for automatic Feynman computation packages is clear as this task is beyond the 
reach of the theorist community. 
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3. Components of an Feynman diagrams automatic computation package 

The general structure of a typical package is represented in the fig.|[ 

3.1. Model definition 

The framework on which the computation takes place has to be precisely and 
uniquely defined. In principle the selection of the Lagrangian should lead to the 
specification of all coupling constants. All fundamental parameters are computed 
from a set of experimental values. Several models have already been implemented, 
including QED, QFD, QCD and MSSM^. For example, QFD can be defined* from 
a selection of the following parameters (a(0), a(Q 2 ), a s (Q 2 ), Mz,Tz, Mw,^w,Gf, 
sin 2 W , Vq KM , Mjj, Mf) depending essentially on the experimental precision. The 
coupling parameters can be computed (Qf,g[, gt,g^,gwwz) as well as M\y and 
rv the W propagator mass and width. Then the final observables (|i?| 2 , d ^ e ,cr tot ) 
should be uniquely defined. 

3.2. Process definition 

The initial and final state particles are selected as well as the order in the coupling 
constant (a,g,a s ). For example: e + e~ — > W + W~~~{ at 0(a 3 ) for tree level or 
e + e~ — » W + W~ at 0(a A ) for 1 loop corrections. The final state can be generic 
like e + e~ — > 4 fermions or 77 — > qq. The initial state may contain composite 
particles as in P7 — » P<y. The way to handle radiative correction (ISR, soft photon 
contribution) and structure functions should also be defined. 

3.3. Graph generation, drawing and selection 

At the tree level, the graph generation is straight forward, the generalization to 
the n-loop case has been solved by the so-called "orderly algorithm"ll2 where all 
possible topologies without duplication are identified. However the computing time 
at a given order 0(a n ) is proportional to the #nodes!§. Performance improvement 
for complex event has been achievedta based on the use of vertex classification. 

The graph drawing module, producing usually postscript diagram representa- 
tion, is needed for visual check, interactive selection and inclusion in publication. 
This work of art has been completed in some case up to 2-loop diagrams where 
topology difficulties becomes quite substantial. 

Diagram ordering or classification based on some physics interest is a very im- 
portant issue. Finding rules based on topological informations and model properties 
to select gauge invariant graph subgroups or to order diagrams by their expected 
contributions would provide ways to reduce computation loads. 

T Quantum Electro Dynamic, Quantum Electroweak Dynamic, Quantum Chromo Dynamic and 
Minimal Supersymmetric Standard Model 
tLEP-II workshop proposal (see R. Kleiss's talk) 
§A node is an external particle or a vertex 
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3.4. Matrix element elaboration 

Having generated all the diagrams and using the model description, one can now 
write the matrix element function. It will provide the contribution of all selected 
diagrams for each phase space point. Two approaches are possible and have been 
used. Let us consider, for example, the inclusive reaction, e + e~ — > X. The matrix 
element can be written as follows: 

|i*l 2 = £lE^'OM40l 2 (i) 

where are helicities and g graphs. 
The symbolic approach 

i^i 2 = E E OM*, o<M4 or w 

g,g' 

\R\ 2 = ^Tr(ff0 g p&g<) where = 7V7 (3) 
g,g' 

is based on trace calculus using symbolic manipulation languages: Schoonschip, 
Reduce, Form, Maple, Mathematica. The sum of the trace for all graphs is per- 
formed after summing over the helicities of the external particles. This technique 
leads to compact expressions for simple processes and to the analytic cancelation 
of singularities. However large expressions are difficult to reduce (nb. terms oc (nb. 
of diagrams) 2 ), polarization effect cannot easily be computed as it adds many more 
terms and the gauge invariant check is difficult as the computation can only been 
done in the most appropriate gauge scheme. 

In the numerical approach, expression ([j]) is directly developed in wave func- 
tions, vertices and propagators components. The matrix element is then built from 
sequential calls to a helicity amplitude library containing all necessary numerical 
routines associated to each basic components. This is a systematic approach, valid 
for all tree level processes and not limited by the complexity of the diagrams. How- 
ever, numerical stability hampers the practical use of this technique. The load is put 
on the integration package and its ability to deal with singularities. Large comput- 
ing time is needed for the integration and event generation. For example 100 hours 
on an HP-735 were necessary to compute a (2 —>■ 5) process like e + e~ — > e~i7 e ud'y. 

3.4.1. Higher order computation 

Three cases should be considered: 1-loop, 2-loop and n-loop corrections. In the 
first mixed approach (symbolic and numerical) is followed. Tree level, 1- 

loop diagrams and counterterms are generated, then the product of loop and tree 
diagrams is performed symbolically: (ioop), * (tree)*. The output is composed 
of 2-4 point functions. A numerical integration is then performed using numericaHj 
or analyticalHl scalar loop integral libraries( up-to 5-point functions). Although 
limited to 4 external particles in QFD, automatic computation have been used for 



8 Towards a Complete Feynman Diagrams Automatic Computation System 



e + e~ — ► tt (2 tree and 50 1-loop diagrams) or e + e~ — » Z°H (1 tree and 80 1-loop 
diagrams). Numerical instabilities and computer performance are the major limita- 
tions of this approach. Two-loop calculations is a very active field where hand com- 
putation are still necessary to find approaches suitable for automatization. These 
calculations imply dealing with quite complex renormalization and regularization 
problems. However, many tools are being developed, many partial results have been 
obtained and systematic approaches begin to emerge. Higher loop calculations are 
still in their infancy, no automatization have been yet tried, only computations with 
the help of algebraic tools has been achieved. Complete automation are probably 
impossible^ 

3.5. Integration and event generation 

The matrix element function is integrated over the multi-dimensional phase space 
restricted by the cuts introduced by the experimental acceptance and by the in- 
trinsic constraints needed to tame the gauge violation divergences. The mapping 
of the integration parameters to the physics variable is called the kinematics. This 
transformation must regularize or (at least) decorrelate singularities (infra-red di- 
vergences, mass singularity, 7-t channel singularities, resonance formation) to mini- 
mize the variance in order to obtain the best numerical stability. Each singularity is 
mapped to an independent variable whenever possible. Otherwise, the phase space 
is split to adapt to each kinematical regions. 

Multi-dimensional adaptive integration algorithms have been developed. In the 
"stratified sampling" approach, the gricLspacing is adapted to the integrand mag- 
nitude and gradient (VEGASlHa, BASESta). Iterative integration is then performed 
until the requested accuracy is reached. 

New ideas based on quadtree and simplexes partitioningS or on waveletS anal- 
ysis are proposed to improve the performance of this major module. 

The gauge invariance is checked numerically by selecting a point in the phase 
space and computing the matrix element for the various gauge parameters and 
gauge schemes (covariant, unitary, axial gauge). 

Structure functions F(x, s) can be taken into account to the cost of increasing the 
number of dimensions in the integration er(s) = J dxF(x, s)<ro{x, s). The kinematics 
becomes also more complex due to the s-dependence of the integrand. 

Parton shower and hadronization are performed by independent packages. The 
partonic final state is transformed into physical particles using one of the many 
hadronization scheme including: Independent fragmentation (COJETS^), Color 
string (JETSET/PYTHIAll, ARIADNE0), low-mass cluster (HERWIgI!). Color 
correlation and helicity constraint (resonance formation) are more difficult issues 
and should be implemented carefully. 

Event generation is obtained by sampling the distributions computed during the 
integration step. SPRINGER is an event generator which has been developed to use 



tSee D. Broadhurst's paper in these proceedings 
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Fig. 3. A Grc++ example 



the BASES integration informations. 



Fig. 4. Standardization: an Object Oriented 
Approach 



4. Packages 

To date' only two packages, GRACE and COMPHEP, are fully automatic: from 
the process definition to the cross-section final values. 

The GRACE packages have been developed by the Minami Tateya Collaboration! 
GRACE- T0 is a general tree level package which adopts a numerical approach based 
on helicity amplitudes for massive particles using the CHANEL library. The QFD 
and QCD model are built-in. Anomalous couplings and MSSM are almost ready to 
be released. Integration and event generation rely on the BASES/SPRING pack- 
age. GRACE-LI is the one-loop extension to the previous package. The generation 
of n-loop and counter terms is performed, but the rest of the package is limited to 
1-loop and 4 external particles. Algebraic computation is performed on the product 
of the loop and the tree level diagrams using Reduce. Then a numerical integration 
involving scalar loop integrals is carried out. Renormalization is performed auto- 
matically in 4+e dimensions. e + e~ — ► tt and 77 — > tt are some of the processes 
which have been calculated with this package. GRC++ is an interactive package 
based on KUIP embedded in PAW++. It is fully interactive for up to 2— >4 and 
some 2— »5 processes (fig .^). 

The COMPHEP package has been developed by the INP Moscow GroupN It 
is a general tree level system (up to 2 — > 4 process). It is based on an analytical 
approach where the square of the sum of the diagrams arc computed using their own 
symbolic calculation package, Reduce or Mathematica. The complete QFD model is 



t Check on these proceedings for more recent updates 
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treated in unitary and 'tHooft-Feynman gauge. Anomalous coupling and MSSM 
are in preparation. Furthermore the user may define its own model. Structure 
functions are implemented using PDFLIBEl BASES and SPRING are used for 
integration and event generation. The kinematics is generated automatically for 
the most current needs. The package provides an interactive computation for the 

2- ^3 and some 2^4 processes: from process specification to energy dependent total 
cross-sections and various angular distributions. 

MADGRAPH is an automatic diagram generatiorSpackage for tree level pro- 
cesses. It produces a matrix element in term of massless or massive helicity ampli- 
tudes using the HELAS library. The QFD and QCD model are available, MSSM 
will be released soon. BRS gauge invariance check is automatically built by the 
package. Diagram generation is interactive but the kinematics, integration and 
event generation are let to the user. 

FEYNARTS/FEYNCALC from the Wurzburg GroupEl is a Mathematica 
package providing convenient tools for radiative corrections in the Standard model. 
It generates and computes all tree, 1-loop graphs and some 2-loop processes. It can 
handle up to 4 external particles (2^2 at 1-loop and 2-loop photon self energy). 
QFD and QCD are basically built in and MSSM is in preparation. The 1-loop 
e + e _ — > HZ calculation^ has been performed with this package. 

FDC, the Feynman Diagram ComputatiorH package performs automatic calcu- 
lation of 1-loop diagrams using the Wick's theorem. It follows an analytic approach 
based on Reduce, RLisp to produce the matrix element. This package is still in 
development stage. 

5. Computing Aids 

Some tools have been developed for more specific computations and are sometimes 
more advanced that some of the modules of the complete packages. However they 
cannot be used blindly to automatically produce cross-sections and event generators. 

hipEI is a high level Maple or Mathematica functions library to perform sym- 
bolic calculation on trace of product of Dirac Matrix. It provides an easy way 
to elaborate cross-section integrand and decay width functions. It implements the 
Vector Equivalence technique where pair of external fermions are replaced by an 
equivalent four- vector. 

TRACER from the Garchin groupEl is a Mathematica package to perform 
symbolic manipulation and trace operations on string of 7-algebra objects in n- 
dimensions usingThe 't Hooft-Veltman scheme. 

COMPUTE^ from the Raleigh GroupEl is a Maple or Mathematica package 
implementing spinor techniques for exclusive processes in perturbative QCD, for 
example: Nucleon compton scattering (378 basic Feynman Diag.). 

MINCErU is a Form package dedicated to the calculation of massless 1, 2, 

3- loop diagrams of propagator type, for example: Z° — > Hadrons (a,, NNL approx- 
imation) . 

PHYSIC aEI is a Mathematica package for the symbolic calculation of tree level 
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processes. 



6. Standardization 

It was decided during the workshop to set up a standardization agreement in order 
to open these packages to the outside world and to fix the structure of the event 
generators. Working groups have been formed and the PISA-I document will be 



available soon on bttp : / /lapphpO . in2p3 . fr /aihep /aihep .html 



General package structure: A unique definition of the input and output of 
each major module should permit the building of a unified workbench (Fig.^J) where 
modules could be exchanged easily and where the same output could serve as input 
to several other modules. 

General structure of the event generator: The event generator routine 
needs a deeper level of standardization in order to be interfaced to various struc- 
ture functions, radiative corrections and hadronization packages. Moreover the gen- 
eral structure should permit a simple introduction in the large detector simulation 
packages. 



7. Conclusions 

Complete tree level automatic computation systems like GRACE, COMPHEP and 
MADGRAPH, complementary in their approaches, are available today. Improve- 
ments in the code efficiency, in the user interface and in the kinematics libraries 
are been pursued. Based on_these packages, cross-sections and event generation 
databases are being preparedcS E3 for all tree level processes(2 — > 2, 3, 4, 5). Com- 
plete 1-loop QFD computations are close to be released. Although numerical sta- 
bility problems and limitation to 4 external particles may reduce its practical use, 
it clearly demonstrate the feasibility of this programme and, judging from the work 
being dedicated to these issues, it will keep improving and getting more general in 
the coming years. QCD packages are now being actively developed. Next to lead- 
ing order computations are investigated in the framework of supersymmetric and 
string based approaches. Complete automatic 2-loop correction packages are still 
in the science fiction section as they are one order of magnitude higher in technical 
complexity and in the need for computer performances. However a lot of activity 
is devoted to these themes and progress will definitely come. Complete automatic 
n-loop calculation seems out of reach with today technology. The future of au- 
tomatic Feynman diagrams automatic computation depends on a widely accepted 
standardization scheme. Working groups have been created during this workshop 
to produce the first first Standardization Agreement document (PISA-I) very soon. 
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